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We report high-resolution specific heat measurements on the quasi-2D organic superconductor k- 
(BEDT-TTF)2Cu(NCS)2, performed by using a most sensitive ac-modulation technique. The main 
observations are (i) a discontinuity at T c much in excess of what is expected for a weak-coupling BCS 
superconductor and (ii) a quasiparticle contribution to the specific heat with an exponentially weak 
temperature dependence at T <C T c . The latter finding is incompatible with an order parameter 
that vanishes at certain parts of the Fermi surface. The data for T < T c can be well described by a 
strong-coupling extension of the BCS theory - the a-model - similar to what has been recently found 
for the K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br salt [Elsinger et al, Phys. Rev. Lett. 84, 6098 (2000)]. 

PACS numbers: PACS numbers: 74.70.Kn, 74.25. Bt, 65.40.Ba 



The quasi-twodimensional organic charge-transfer salts 
k-(BEDT-TTF)2X based on the electron-donor molecule 
bis(ethylenedithio)-tetrathiafulvalene, commonly abbre- 
viated to BEDT-TTF or simply ET, are of particular 
interest because of their rich phase diagram including 
superconductivity at relatively high temperatures next 
to an antiferromagnetic insulating state |l], 0. The 
pressure- induced transition from an antiferromagnetic in- 
sulator to a superconductor for X = Cu[N(CN) 2 ]Cl is 
one of the intriguing observations that places these ma- 
terials in line with other exotic superconductors where 
pairing mechanisms which are different from the con- 
ventional electron-phonon interaction might be of rele- 
vance. In fact, the proximity of magnetic order and su- 
perconductivity along with the presence of spin fluctua- 
tions above T c as inferred from 13 C-NMR measurements 
H |J, |s| , have been taken as strong indications for a spin- 
fluctuation-mediated superconductivity || [?], g, ||, [R| 
similar to that which has been proposed first for the 
high-T c cuprates |n]]. Since the identification of the 
relevant pairing interaction is a very difficult problem, 
many experiments have focussed on the determination 
of the structure of the superconducting order parame- 
ter - the gap amplitude A(fe) - which is intimately re- 
lated to the pairing mechanism. However, despite in- 
tensive experimental efforts devoted to this issue, no 
consensus has been achieved yet |Q. Arguments in 
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favour of an unconventional order-parameter with d- 
wave symmetry for the above k-(ET)2X salts have been 
derived from temperature-dependent investigations, no- 
tably NMR measurements @, [l3], [13], thermal conduc- 
tivity pl| and one specific heat study E3|. More recent 
attempts such as STM [fl7| , mm- wave transmission or 
thermal conductivity |1S[ have focussed on orientational- 
dependent investigations aiming at a direct determina- 
tion of the gap anisotropy. Although a modulation of 
the gap structure has been seen in all of these studies, 
they come to quite different conclusions on the direction 
of the gap zeroes, see e.g. Q. 

Conversely, there are numerous experimental studies that 
are consistent with a conventional BCS-type of supercon- 
ductivity. Among them are measurements of the mag- 
netic penetration depth [^0) and surface impedance |2^| 
as well as the observation of a BCS-like mass isotope 
effect [^2| and a pronounced superconductivity-induced 
phonon renormalization |25)| . The latter two experi- 
ments clearly demonstrate the relevance of intermolec- 
ular phonons for the pairing interaction. 
A very powerful method to probe certain aspects of the 
gap structure - in particular the question whether gap 
zeroes exist or not - is provided by specific heat measure- 
ments. In case this integral thermodynamic technique 
were to find a low temperature electronic quasiparticle 
contribution, C' cs , that varies exponentially weakly with 
the temperature, the existence of gap zeroes on the Fermi 
surface could be definitely ruled out. On the other hand, 
the observation of a non-exponential temperature depen- 
dence does not represent a clear proof of gap zeroes as 
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it may also originate in extraneous contributions such 
as impurity phases, normal-conducting regions or pair 
breaking. In fact, an exponential low-T specific heat be- 
haviour implying a finite energy gap has been found in 
recent high-resolution specific heat measurements on the 
X = Cu[N(CN) 2 ]Br salt [§4|. Moreover it has been shown 
in this study that the T 2 dependence in C cs reported for 
the same compound by Nakazawa et al. Jl6| most likely 
originates in their incorrect determination of the phonon 
background. 

Here we present high-resolution specific heat measure- 
ments on the related k-(ET) 2 Cu(NCS) 2 salt with the in- 
tention to clarify the existence or absence of gap nodes 
for this compound. 

The single crystals used were synthesised by the stan- 
dard electrocrystallization technique as described else- 
where 



25 



For the present investigations, two high- 
quality single crystals with regular plate-like shapes and 
masses of 0.63 mg (# 1) and 0.72 mg (# 2) have been se- 
lected. The specific heat has been measured utilising a 
high-resolution ac- modulation technique |2(|. The setup 
has been designed for investigating very small plate-like 
single crystals such as the present compounds. The sam- 
ple holder consisting of a resistive thermometer (Cer- 
nox CX-1080-BG) and heater is attached to a 4 He-bath 
cryostat equipped with an 8T superconducting solenoid. 
The calorimeter has been checked by measuring high- 
purity samples of Cu and Ag with typical masses of about 
3 ~ 4mg which have absolute heat capacities compara- 
ble to those of the small crystals studied here. In the 
temperature range 2K < T < 30 K, the maximum devi- 
ations from the literature results amount to ±2 %. 

Fig. ^ shows the specific heat as C/T vs T of crys- 
tal # 1 over the entire temperature range investigated 
p7[ . The phase-transition anomaly at T c is clearly vis- 
ible although it amounts to only about five percent of 
the total specific heat at this temperature. To determine 
the quantity of interest - the quasiparticle specific heat 
in the superconducting state C' cs - one has to get rid of 
the large phonon background C pn . It is known for these 
molecular systems that, due to low-lying optical phonon 
modes, C pn starts to deviate from a Debye-like behaviour 
already at low temperatures 
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ties related to C pn , we proceed in analyzing the difference 
AC(T) = C(T, B = 0)- C(T, B = 8T). For the field ori- 
entation used in our experiment, B \\ a* , where the upper 
critical field is B C2 « 6 T J3TJ] , the data taken at 8 T rep- 
resent the normal-state specific heat C n . Provided that 
C n is field independent, which has been found in previ- 
ous measurements for the present compound |nj and also 
for X = Cu[N(CN) 2 ]Br |j, the quantity AC has the ad- 
vantage that the unknown C p h and all other extraneous 
contributions (such as a small amount of Apiezon grease 
used to improve the thermal contact to the thermome- 
ter and heater) cancel each other out. The experimental 
fact that C n is independent of the magnetic field within 
the experimental resolution |24|, |3lj is consistent with 
C n = C p h + 7? 1 p2l , where 7 is the Sommerfeld coeffi- 



cient. 

Fig. H shows the specific heat data for B = and 
B = 8 T applied along the a*-axis, i.e. perpendicu- 
lar to the conducting planes. For technical reasons, 
measurements in fields of B = 8 T were limited to 
T < 8K. For C{T,B = 8T) = C„(T) at T > 8K, 
an interpolation based on a polynomial fit between the 
C(T < 8K,B = 8T) and C(T > 10K,B = 0) data 
was used, cf. the solid line in Fig. |^. The so-derived 
AC = C(0T) - C(8T) = C os - lT of crystal # 2 is 
shown in Fig. ^ together with AC expected from the 
BCS weak-coupling theory |33| . The theoretical curve is 
based on a Sommerfeld coefficient 7 = (23±1) mJ/molK 2 
as determined by low-temperature C(T) measurements 
employing a thermal-relaxation technique on a compila- 
tion of three single cry stals with total mass of 25.8 mg, 
see inset of Fig. [y |34|j . This 7 value is consistent with 
7 = (25 ± 3)mJ/molK 2 derived from an earlier specific 
heat study on the same compound f3l|| by analysing the 
data in the range 1.3K<T<3K which is slightly above 
that used here fl35| . 

Fig. ^ demonstrates that, similarly to what has been 
recently observed for the X = Cu[N(CN) 2 ]Br salt ||, 
AC(T) deviates markedly from the weak-coupling BCS- 
behaviour in both the jump height at T c as well as the 
overall temperature dependence. However, as for X — 
Cu[N(CN)2]Br, a much better description of the data 
is obtained by using the semi-empirical extension of the 
BCS formalism to strong-coupling superconductors - the 
so-called a-model |pq| . It contains a single free parame- 
ter a = A(0) /ksT c [ks being the Boltzmann constant) 
which scales the BCS energy gap A(T) = (a/aecs) ■ 
Abcs(^ 1 ) with aecs = 1-764. As Fig. || clearly demon- 
strates, the strong-coupling BCS model with a = 2.8±0.1 
provides an excellent description of the data over the en- 
tire temperature range investigated. The error margins 
in a account for the uncertainties in both 7 as well as 
T c . In the inset of Fig. || we show AC data of crystal # 1 
together with the theoretical results of the a model. Al- 
though T c of crystal # 1 is slightly reduced compared to 
that of crystal # 2, we find again an excellent agreement 
with the strong-coupling results employing the same pa- 
rameter a — 2.8. 

The fact that both data sets in Fig. ^ are well de- 
scribed within the strong-coupling BCS model implies 
(i) that C es reveals an exponentially weak temperature 
dependence at low temperatures and (ii) the thermody- 
namic consistency of the data, i.e. entropy conservation. 
The exponential variation of C es becomes clearer in Fig. 
^, where C cs /^fT c is shown in a semi-logarithmic plot 
as a function of T c /T. Here, C cs = AC + 7T with 
7 = (23 ± l)mJ/molK 2 has been used. The solid line 
represents the same strong-coupling curve shown in Fig. 
||. The figure also includes the weak-coupling BCS result 
(dashed line) which predicts C es /jT c oc exp (—a/^T c /T) 
with a A = 1-44 for 2.5 < T c /T < 6 |7). The data 
of k-(ET) 2 Cu(NCS) 2 presented here follow C cs /jT c oc 
exp(— 2.5T C /T) down to the lowest accessible tempera- 
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ture. The enhanced prefactor in the exponent of a& = 2.5 
as compared to 1.44 for the weak-coupling BCS model 
reflects the strong-coupling character of the present su- 
perconductor. 

Fig. [| demonstrates that our results are fully consistent 
with an exponentially vanishing C cs at low temperatures 
and, thus, an energy gap without zeroes at the Fermi sur- 
face. The same conclusions have been drawn from a sim- 
ilar analysis of C(T, B) data on K-(ET) 2 Cu[N(CN) 2 ]Br 
p4j . As the specific heat is an integral technique that 
picks up all excitations at the Fermi surface, the above 
findings are incompatible with the existence of gap ze- 
roes as claimed by other experiments, notably NMR 
H |l3| . Since the latter results have been obtained in 
a finite magnetic field applied parallel to the conducting 
planes, the discrepancies are possibly related to an as yet 



not understood influence of the magnetic field. 
In conclusion, high-resolution specific heat measure- 
ments on small high-quality single crystals of k-(BEDT- 
TTF)2Cu(NCS)2 have been performed in both the nor- 
mal and superconducting state. The data analysis em- 
ployed, which minimises uncertainties associated with 
the unknown, large phonon background reveals an elec- 
tronic quasiparticle contribution to the specific heat that 
varies exponentially weakly with temperature at T <C T c . 
This behaviour is fully consistent with a finite energy 
gap all over the Fermi surface and rules out the exis- 
tence of gap nodes. Moreover, we find an excellent agree- 
ment with the predictions of the strong-coupling variant 
of the BCS model employing an a parameter 2.8 that 
slightly exceeds the one found recently for the related 
K-(ET) 2 Cu[N(CN) 2 ]Br salt ||. 
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FIG. 1: Specific heat data as C/T vs T of fi>(ET) 2 Cu(NCS) 2 , 
crystal # 1. The inset shows data as C/T vs T 2 on a compi- 
lation of several single crystals with a total mass of 25.8 mg 
taken at low temperatures in a magnetic field B — 10 T. The 
solid line is a linear fit of the form C/T = J + /3T 2 to the data 
below T = 2 K. 




FIG. 2: Specific heat data as C/T vs T of k>(ET) 2 Cu(NCS) 2 , 
crystal # 2, in the vicinity of the superconducting transition. 
Data have been taken in zero field and in an overcritical field 
B = 8 T. The line represents a polynomial fit to the data 
C(T < 8K,B = 8T) and C(T > WK,B = 0), see text. 




FIG. 3: Specific heat difference AC = C*(OT) - C(8T) = 
C(OT) - Cn of crystals # 2 and # 1 (inset). The dotted and 
solid thick lines represent the BCS curves for weak and strong 
coupling, respectively. 




FIG. 4: Semi-logarithmic plot of the electronic contribution 
to the specific heat in the superconducting state as C cs /"/T c 
vs T c /T. The dotted and solid lines represent the weak- and 
strong-coupling BCS behaviour, respectively, see text. 



